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Abstract 

The reconstruction approach [Shu C.W.: SIAM Rev. 51 (2009) 82—126] for the numerical approxi- 
mation of f'(x) is based on the construction of a dual function h(x) whose sliding averages over the 
interval [x - ^Ax, x + ^Ax] are equal to /(x) (assuming an homogeneous grid of cell-size Ax). We study 
the deconvolution problem [Harten A., Engquist B., Osher S., Chakravarthy S.R.: J. Comp. Phys. 71 
(1987) 231-303] which relates the Taylor polynomials of h{x) and /(x), and obtain its explicit solution, 
by introducing rational numbers t„ defined by a recurrence relation, or determined by their gener- 
ating function, grix), related with the reconstruction pair of e*. We then apply these results to the 
specific case of Lagrange-interpolation-based polynomial reconstruction, and determine explicitly the 
approximation error of the Lagrange reconstructing polynomial (whose sliding averages are equal to 
the Lagrange interpolating polynomial) on an arbitrary stencil defined on a homogeneous grid. 
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differences, finite volumes 
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1. Introduction 

The Godunov approach lU to hyperbolic conservation laws 

d,u + d^Fiu) = (1) 

is based on space-time averaging of the PDE ([l). Assuming an homogeneous time-independent grid 
(Ax = const), space-averaging of ([1]), over the interval [x - ^Ax, x + |Ax], leads to the exact relation yj] 



5 

dt"'^"''' ' Ax 

where 



u(x, t) + ^[F (u(x + jAx, O) - F (u(x - jAx, f))] = (2) 
u(x,t):- I u(x + ^Ax,t)d^ (3) 



are the sliding cell-averages of the solution. Defining the sliding cell-averages F(u), by applying the 
operator ([3) on F(u), we have immediately by differentiation, provided that Ax = const, 

dF(u(x,t)) _ + i Ax, f)) - f[u(x - iAx, t)) 

dx Ax 
exactly, so that, combining I© and @ 

d,u + d^Jiu) = (5) 
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ie the equation for the sHding cell-averages, for Ax = const, has the same form as the original equa- 
tion lH. For this reason, it is assumed that what is computed (and stored at the nodes of the computa- 
tional grid [3, 4]) are the cell-averages of the solution. 

In the discretization of© we are led to consider the computation of the derivative of a function f(x) 
(corresponding to S) sampled on the computational grid, by differences at x + ^Ax of the values of an 
unknown function h(x) (corresponding to m), which has to be reconstructed iSSSHBl from the values 
of its cell-averages sampled on the grid. In the following, we concentrate on the spatial discretization 
problem, viz compute f'(x) via reconstruction ofh(x ± ^Ax) B 

Reconstruction (Definition Ol is the basis of ENO 0, E E Ml and WENO 10, H, H, H O, [H, [H, 
13, fisll schemes. Although the term polynomial reconstruction is quite general, including Hermite- 
interpolation |16], spectral methods [17], and spectral element techniques fl8^, we concentrate in 
the present work on reconstruction approaches based on Lagrangian interpolation of the function 
averages. There exist several algorithms for Lagrangian-interpolation-based polynomial reconstruc- 
tion and these have been successfully used for the construction of progressively higher-order 
schemes i6[[i3,[l^, using symbolic calculation |13, 14]. The reconstruction via primitive approach |0] 
is probably the most widely used in order-of-accuracy proofs lIMSl, while the reconstruction via decon- 
volution approach [2] has been formulated with respect to the solution of linear systems. Most of these 
schemes and associated order-of-accuracy relations LS, XL 2A, iSJ were developed for particular 
values of the order-parameter r (determining the discretization stencil), using symbolic computation. 
On the other hand, analytical relations for the order-of-accuracy of the approximation of h{x), for arbi- 
trary reconstruction-order-parameter r are not available. To obtain such relations it seems necessary 
to study in detail the relations between a function h{x) (which is reconstructed) and it^ cell-averages 
/(x). This is the reconstruction via deconvolution approach defined by Harten et al. |2]. Up to now, 
these relations were obtained by solving, using symbolic calculation, the associated linear system 10, 
(3.13b), p. 244], up to a certain order. Although the solution by symbolic computation of the linear 
system |2, (3.13b), p. 244] is not difficult, it is only valid up to a certain 0(Ax'^), and the non availability 
of an explicit solution hinders the development of general expressions of the approximation error of 
the reconstruction. 

Along with the numerous successful developments of practical WENO schemes based on the re- 
construction via primitive approach 10, Q], the unknown function h{x) which is reconstructed by its 
cell-averages /(x) appears explicitly in recent analyses ifisl. [l3] of the truncation error. Analyzing 
the reconstruction error in terms of the unknown function h{x) and its derivatives {reconstruction via 
deconvolution [2]) is a more intuitive approach, especially when considering the discretization error 
of the WENO ^_approximation to fix) and potential improvements in the formulation of the nonlinear 
weights fis", 14]. Analyses based on the reconstruction via primitive |2] approach are somehow less 
straightforward as they involve the primitive of the reconstructed function £^ h{()d^^ The main mo- 
tivation of the present work is to contribute to the development of analytical tools, applicable to the 
study of the truncation error flsL [l4ll . determination of the loss-of-accuracy at smooth extrema |13] 
and research for the improvement of WENO schemes |14], maintaining the in-built scalability (with 
the stencil width) towards higher order-of-accuracy of WENO schemes lla JlUllSn . For these reasons, in 
the present work we are not interested in the development of a new algorithm for the solution of the 
reconstruction problem. Instead, we focus on reconstruction relations of general validity, ie stencil- 
independent, which are necessary for the study of the approximate reconstruction order-of-accuracy. 

In !j2]we study the general relations underlying the reconstruction approach for the numerical ap- 
proximation of the 1-derivative /'(x) of a function /(x). Initially we study the relations between the 
derivatives of a function /(x) and those of a dual function h{x), whose sliding averages, over a con- 
stant length Ax, are equal to fix). We will call the functions, /(x) and /i(x), satisfying this relation 
a reconstruction pair for the discretization of fix) (Definition 12. 11 1. We introduce the rational num- 
bers T„ € Q, defined either by a recurrence relation (Lemma I2.5I > or through a generating function 



^The lower bound of the integral being of no particular consequence I3ll3. ll9r we can chose the coordinate of grid-point xi 
instead of the usual (but more abstract) lower bound at -oo Ia.l4. ll9r . 
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(Theorem l2.9l l, which are used to develop exphcit series representations of h(x) (and of its derivatives) 
with respect to powers of Ax and the derivatives of f{x). The principal new result in Sj2] is that we 
are able to give explicit solutions to the fundamental relations of the reconstruction via deconvolution 
approach ^ (3.13), pp. 244-246], which (Lemma l2.5D are widely used throughout the paper. The gen- 
erating function of the rational numbers t„ e Q appears in the expression of the reconstruction pair of 
(Theorem[221l. 

In ^ we study the particular case of polynomial reconstruction. We show (Lemma 13. ID that for 
every polynomial p j (x) of degree M in x we can define, using the numbers t„ (Lemma l2.5l l, a polynomial 
Ph(x), also of degree M in x, so that pf and pi, are a unique reconstruction pair (Definition l2.1D . Initially 
( ^3.21 1 the numbers t„ (Lemma 12.51 1 were introduced, using a matrix algebra approach to study the 
relation between pf(x) and pi,(x)- This part of the paper ( ^3.21 1 gives the explicit inversion of the matrix 
appearing in the reconstruction via deconvolution theory [2, (3.13b), p. 244]. 

In practice f(x) is usually approximated by its Lagrange interpolating polynomial pf(x; S,,m_,m+, Ax) 
on a given stencil S,,m_,m+ (Definition 14.11 1. and /i(x) is approximated by the reconstruction pair of 
Pf(x; S,;m-,m+, Ax), ph{x; S,;m_,m+, Ax), which we will call the Lagrange reconstructing polynomial fM- In 
fjllwe study the approximation error of the Lagrange reconstructing polynomial, Ei,(x; S,,m_,m+, Ax) :- 
Ph(x; S,,M_,M+, Ax) - h(x), and obtain an explicit relation for the expansion of this error in powers of Ax 
(Proposition 14. 7D . This is only possible through the explicit solution of the deconvolution problem 
(Lemma 12. 5D . In fj5] we briefly summarize the existence and uniqueness results concerning the re- 
constructing polynomial. Finally, in fj6]we briefly describe some applications of the present results to 
practical WENO schemes, and discretization methods in general, highlighting the merits of the recon- 
struction via deconvolution approach, as developed in the present work. 

Standard results referring to the Lagrange interpolating polynomial ||20| . l2l}] are included only 
when they are necessary for the proof of the new results concerning the reconstructing polynomial. 
Useful relations for summation indices in multiple sums ll2^l23ll . and other identities, used throughout 



the paper, are summarized in Appendix A 



2. Reconstruction pairs and exact reconstruction relations 

Before proceeding to a detailed examination of the reconstruction of polynomials we examine the 
general relations underlying the reconstruction approach for the evaluation of the derivative /'(x) of a 
function /(x), via the construction of a function h(x) (reconstruction pair of /(x); Definition l2.1l l. whose 
sliding (with x) averages |3, 4] on the interval [x - jAx,x + jAx] are equal to /(x), over an appropriate 
interval x e / c R. We express in particular the derivatives of h(x) as series of the derivatives of 
fix), with coefficients determined by the derivatives at Ax = of the function gAAx) appearing in the 
reconstruction pair of the exponential function (Theorem l2.9l >. 

2.1. Reconstruction pairs 

The basic idea underlying reconstruction procedures to compute the derivative /'(x) of a function 
fix) follows directly from the Leibniz rule [24, pp. 411-412] giving the derivative of a definite integral 
with respect to its (variable) bounds. To this end we need to construct a function h(x) whose sliding 
(with x) average over an interval [x - ^Ax, x + jAx] of constant width Ax is equal to f(x). 

Definition 2.1 (Reconstruction pair). Assume that Ax e R>o is a constant length, and that the 
functions / : / — > M and h : I — > M are defined on the interval / = [a - jAx, b + jAx] c M, satisfying 
everywhere 

/(x) = — h(Od( \/xe[a,b] (6a) 

Ax J.,_^Ax 

assuming the existence of the integral in ( I6al l. We will note the functions /(x) and h(x) related by i6ai 

h =^(l;A.v)(/) (6b) 

/ (60 
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and will call / and h a reconstruction pair on [a, b], in view of the computation of the 1-derivative. □ 
By definition dHal l. R^^.^^y ( l6cl ) is defined by 

Vx e [a, b] 

AJC Jx-Ux 



'^x+jAx 



( f^x+'-Ax ] (7) 

\h:[a-\Ax,b+\Ax] — h{Od( eR'^xe[a,b]\ 



ie ^(I'^j.) is a mapping applicable to all real functions defined on [a - jAx,b + jAx] c M for which the 
integral ( l6al ) exists. Then, R(i-ax} is defined as the inverse mapping of R^^^.^^,^ assuming that the 
inverse mapping exists (cf Remark l2.8l >. 

Lemma 2.2 (Reconstruction). Consider the functions f(x) and h(x) constituting a reconstruction pair 
on [a,b] c E (Definition l2.lt . Assume that f(x) and h(x) are of class (N e N) on the interval I - 
[a - ^Ax, b + jAx] c R. Then 

h^"-^'>(x + iAx) - h^"-^\x - iAx) 

f"^(x) = ^ 2 \/x e [a,b] Vn e {1, ■ • • , A?) (8) 

Ax 

Proof. Direct differentiation of ( l6aD . yields 

h(x + iAx) - h(x - iAx) 

/'(x) = ^- — —— ? — - Vx € [a, b] (9) 

Ax 

by application of the Leibniz rule lHJ, pp. 411-412], and taking into account that Ax is constant Vx. 
Successive differentiation of (HJl yields □ 

All reconstruction-based approaches |0,l2,[l,[l,[l3,IH,[ll,[2l,[l3] for the numerical approximation 
of PDEs are based on, or can be shown to be related to. Lemma [2. 2 1 These relations ^ are exact rela- 
tions concerning the continuous functions / and h. When /(x) and h{x) are numerically approximated 
consistently, ie in a way satisfying (HJl up to a given order Ax*^+' , then © are satisfied up to some order 
<M+ 1. 

Definition 2.3 (Lagrange reconstructing polynomial). Let /?/ be the Lagrange interpolating poly- 
nomial of the function / on the arbitrary stencil {/ - M_, ■ • ■ , / + M+) of M + 1 equidistant points (M :- 
M_ + M+) around point /. Its reconstruction pair (Definition 12.11 * will be called the Lagrange recon- 
structing polynomial on the stencil {/ - M_, ■ ■ • , / + M+}. □ 

Remark 2.4 (Homogeneous grid). The basic relations underlying reconstruction, which are given 
in Lemma r2.2[ hold iff Ax = const, ie, when used as basis for the numerical approximation of /'(x), these 
relations are only applicable on a homogeneous grid. In the case of an inhomogeneous grid, where the 
spacing Ax(x) is a function of position (Ax : R — > M>()) these relations should be modified to include Ax' 
and ((9a.v/!)Ax'. The general case of an inhomogeneous grid requires specific study. □ 

2.2. Deconvolution 

Obviously, the relations between / and h (Lemma 12. 2D imply that the Taylor-polynomials of /(x) 
can be expressed with respect to the derivatives h^"\x ± jAx), which can themselves be replaced by 
Taylor-polynomials of /z(x). We have 

Lemma 2.5 (Deconvolution oth- /?(1;Aa)(/))- Let fix) and h(x) = [R(UAx)(f)](x) be a reconstruction 
pair (Definition 1211), satisfying the conditions of Lemma I2l2] Then VA^tj e N : A^tj < - 1 



f=0 



/ W- 2^ 22^(2^+1)!'' (^) + 0(Ax ^ v„6No:«<A^-2L^J 



Inversely, 



e=o 



Vx e [a,b] 
Vn eNo : n < A^-2L^J 



where the numbers th (Tab. [T) are defined by the recurrence relations 



k-l k 

^" ~ ' ^^^^ 22*-2' (2;t - 2s + 1)! 2^" (2s + 1)! 



;t>0 



(10b) 



(10c) 



Proof. Approximating h({) (which was assumed to be of class in Lemma [2.21 1 in ^ by the corre- 
sponding Taylor-polynomial (Taylor-jet) of order A^tj 1.27. pp. 219-232] around ( - x yields, VA^tj g N : 
Nt,<N-1, 



Ax ^! 



1 rlAx^f ^1 



J_^/_Ax^+i 1- (-1)^+1 
Ax 2j 9f 



Ax^\2''(<^+1)! 2 



+ 0((^-x)'^"+') 
t/f + 0(Ax'^"+') 
/z<^'(x) + 6>(Ax^"+') 

h^'^\x) + 0(Ax^"+^) Vx e [a, ^] 



and since e No 



we obtain 



e+l^2k+l (/teNo) ^1 -(-1)^+1 =2 



^+1=2A: (/teNo) ^l-(-l)^+i=0 



I %J I 

^ ^ ^ Ax^^ 



-M2f)(x) + 0(Ax2L^J+2^ 



(11) 



(12a) 
(12b) 



(13) 



22^(2^+1)!- 

which is OlOal * for « = 0. Successive differentiation of 0131 1 by x yields OlOaD . 

To invert dlOaP we search for numbers T2s is e No) satisfying VMtj e N ; Mtj < - 1 and Vn e No : n < 

2 

Mtj 

/i(«)(x) =2t2,Ax27(«+2^'(x) + 0(Ax2'^™+2) 



A?-2L^J 



1=0 



y y — /,(«+2.+20(_^) + o(Ax2^-+2) 
y y 42£^^ — M'-2-20(^) + o(Ax2^"+2) 

ZiZj\22^(2^+l)! ' ') ' 



T2,A^' + 0(Ax2'^"+2) 



2Mtj ( min(A:,MTj) 

= Z E 

\^-max(0,^:-MTj, 

=zz 

k=0 \ .5=0 



T2s 



22i-2. (2k-2s+\)\ 



Ax^" /i("+2*Xx) + C(Ax2'^^'+0 



T2s 



22*-2.5 (2A:-2s+ 1)! 



Mtj+2x 



(14) 



because of iA.Si . ( fl4l ) holds, provided that {6ko is the Kronecker 6) 



=0 



which is satisfied if the numbers T2k are defined by dlOcH . Truncating ( I13I I to 0(Ax^^ "J) yields dlObl ). 
The remainder in dlOal l and dlObl ) is 0(Ax^^~2^^^~), because only even powers of Ax appear in these 
expressions. □ 



Table 1: Numbers t„ )19c) satisfying recurrence flOct, for0<n<21. 



To= 1 
Ti = 



T2 



24 

T3 = 



7 

5,760 



T4 = 

T5 = 

-31 

^6 - 967,680 
T7 = 

127 



''"8 ~ 154,828,800 
T9 = 
-73 



^10 ~ 3,503,554,560 
Til = 



1,414,477 



''"12 - 2,678,117,105,664,000 
T13 = 

- -8,191 
"14 - 612,141,052,723,200 

T15 = 

16,931,177 



49,950,709,902,213,120,000 
T17 = 

-5,749,691,557 



''"18 ~ 669,659.197,233,029,971,968,000 

T19 = 

_ 91,546,277,357 

"20 - 420,928,638,260,761,696,665,600,000 

T21 = 



Remark 2.6 (Relation to previous work 28]). The results in Lemma expressing the deriva- 
tives of the sliding cell-averages f(x) with respect to the derivatives of the function h{x) - [R(\ ;Ax){f)]{x), 
are straightforward. In particular dlOat corresponds to JH, (2.15), p. 299]. The new results of 
Lemma 12.51 are the inversion relations dlObl l, which are based on the introduction of the numbers 
T„ dlOcD . These results are the general explicit solution of the linear system written in Harten et al. @, 
(3.13b), p. 244], and provide the exact deconvolution relation between f(x) and [R(i;&x)(f)](^) (Defini- 
tion [ITl), in the case of a homogeneous (A.ji: - const) grid. The general case of an inhomogeneous grid 
requires specific study. The inversion relations dlObl ) are the main building block of the present work, 
as far as error analysis of the reconstruction is concerned. We will show that the numbers t„ dlOcl l can 
also be defined by a generating function (Theorem l2.9l l. □ 
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Corollary 2.7 (Taylor-polynomial of h(x + ^Ax)). Assume the conditions of Lemma [23] Then 



LfJ 



h(x + ^Ax)^J] ZSr' 

,5=0 [i=o 



■s-lf 



f^(.-2^)!j 



Ax' f '\x) + (9(Ax^"+>) 



Proof. Since 



+ 2 = 



n + 1 
n + 2 



Vn = 2fc - 1 



(16) 



(17a) 



jlObl ) can be rewritten as 

Ax"'h^"'\x) 



f=0 



{m + 2€)\ 



(17b) 



In that form ( Il7bl > we have a relation between the coefficients of the Taylor-polynomials of f{x + ^Ax) 
and oihix + ^Ax), expressed in powers of ^. In particular, using ( Il7bl >, we have 



/j(x + ^Ajc) = ^ 



Ax'" h^"'\x) 



(9(Ax^™+') 



A'tj J 
m=0 f=0 



T2f Ax'"+2^ /('"+20(x) 



r + C(Ax^"+i) 



i=0 f=0 



{s-2€)\ 



where we used ( IA.3I > and ( IA.2I ). and the fact that A^tj + 1 < 2L^J + 2. This completes the proof 



(18a) 

(18b) 

□ 



This expression ( I16D is useful in computing the error of numerical approximations to h(x) (Proposi- 
tionSDl. 

Remark 2.8 (Existence and uniqueness). From Definition l2.1l it follows immediately (proof by con- 
tradiction) that every reconstruction pair h - R(i-Ax)if), with h(x) continuous, if it exists, is unique. For 
every h(x) analytic in / with radius of convergence rc^(x), the series dlOaH with n - converges, as 
Ntj — » oo, VAx 6 (o, 2rc,,(-^)), so that (because of uniqueness), for every analjf^ic function h(x) there 
exists a unique function / - R^i.^^^ih). Whether the converse is always true, is an open question. 
Assuming /(x) analytic in / with radius of convergence rc/x), does not automatically imply the con- 
vergence of dlObl ) with n = as A^tj — > ^, because lim„^oo(T2„(2n)!) = oo. The necessary conditions of 
existence require further study. Nonetheless, since lim„^oo t2„ - (Tab. [D and t2„t2„+2 < V« e No 
(Tab. [Hi, the class of functions /(x) for which dlObl l with « = is convergent as A^tj — » °° is not empty. 
It is easy to verify that most of the basic functions /(x) have reconstruction pairs h - /?(1;Aa)(/), as 
do all polynomials of finite degree ( ^3. ID . Whenever any of the series ( [TOl l converges as A^tj — > 
the upper limit of the sums can be readily replaced by oo, to yield complete converging expansions 
(power-series). The Godunov approach [1] to hyperbolic conservation laws d,u + d^F - ([D, is based 
on space-time averaging of the PDE ([1]), to obtain the corresponding PDE, d,u + dxF - ©, for the 
cell-averages li ©. Therefore, with respect to the notation used in Definition 12. 11 u corresponds to / 
and u corresponds to h. In the context of reconstruction procedures ll2,H,SSr4] for the discretization 
of hyperbolic conservation laws, the existence of the solution (integrable function) u (ie h) is assumed, 
so that the existence of the sliding-averages S (ie /) follows (Remark l2.8D . Hence, the results obtained 
in ij2] (where the existence of h is assumed) are directly applicable to the Godunov approach for the 
numerical computation of hyperbolic conservation laws. □ 



7 



2.3. Generating function ofT„ and the reconstruction pair ofe\p(x) :- e " 

As mentioned above (Remark 12.81 * most of the basic functions have reconstruction pairs. The re- 
construction pair of the exponential function plays an important role in the reconstruction relations 
(Lemma l2.5D . because it defines the generating function of the numbers t„ (Tab. [Hi- 
Theorem 2.9 (7;(i;A.t)(exp)). The reconstruction pair ofexp(x) :- e-* is 

[7?(i;A.)(exp)](^) = . ' . = §.(Ax)e^ (19a) 
sinh J Ax 



where the function 



8r(x) := (19b) 
sinh jx 



is the generating function of the numbers t„ (Tab.[T]l satisfying dlOcl l 



Furthermore 



Tn V« E No (19c) 



T2„+i 7T^7TT5?"^''(0) = Vn e No (19d) 
(2n + 1)! 



Proof. From dlObl t. since e ' is of class C°°, we have VA^tj e N 

WtJ ,2« ( ^TJ 

[7?(i;A,)(exp)](x) = 2] ^2"^^'"^^' + = Z 



Vn=0 



e-' + 0(Ax2'^"+2) (20) 



Since lim„^co t2„ = and t2„t2„+2 < V« e No, the alternating {Ax^" > V« e No) series in ( l20l l converges 
as A^Tj — > at least VAx e (0, 1). Defining the function gAx) 



«=0 

suggests that 3 : R — > K such that 

[^(i;A.v)(exp)](x) = g,(Ax)e" (22) 



Using (ISll) in ( [6al 



gives 



= 1 r''%,(Ax)ef^/^ = -^^r(Ax) (e-^^^-' - e-^*^) 
Ax J.,_i Ax Ax ^ ^ 



(23) 



Ax 7 Ax 

g.(Ax) = — — = (24) 

eiA'' - e-i*-^ sinh ^Ax 

proving ( ll9aD . It is a simple exercise to show that the function grix) ( Il9bl l is continuous at x = 0, and 
has continuous derivatives of arbitrary order at x = 0, satisfying 

griO) =1 (25a) 
^P«+i)(0) =0 V« e No (25b) 

Comparing the Taylor-series of grix) ( I19bl > with the series definition of grix) ( l2l] >. and taking into ac- 
count ([25) proves ( ll9cD . ( I25bl > yields ( I19dl >. □ 
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3. Reconstruction of polynomials 

Reconstruction of polynomials (Definitionl2.lt is the basis of ENO |0,[2l and WENO llsl [TH [H. 
14. Iisn reconstructions. We investigate in detail the coefficients of polynomial ( il3.1l l reconstruction 



(26a) 



pairs (Definition l2.1D . 

3.1. Polynomial reconstruction pair 

In this section we consider the case where either f(x) or h(x) in Definition [211] is a polynomial. 

Lemma 3.1 (Polynomial reconstruction pair). Let ph(x, x„ Ax) be a polynomial of degree M 

Ph{x; Xi, Ax) c/,,„ {—^^ 

Then pf{x; Xi,Ax) defined by (Definition 12 .11 1 

p f(x; Xj, Ax) : = — I phU; Xj, Ajc)dl (26b) 

is a polynomial also of degree M, with coefficients cy;„ which can be computed from the coefficients c/,^^^ of 
Ph(x; Xi, Ax) 



/7/(x;x,-,Ax) = £]c/„(^^^) 



m=0 

k=0 



22*(2fe+l)\ 2k 

(m + 2k)\ 
22* (2/t+l)!^''"'*^' 



Vm E {0, • ■ ■ , M} 



Vm 6 {0, ■ ■ ■ , M} 



(26c) 
(26d) 
(26e) 



Inversely, the coefficients ci,^^^ of ph(x; Xi, Ax) can be computed from the coefficients cj^^ of pf(x;xi,Ax) 



1 

y c (m + 2^) ! Vm e {0, ■ ■ ■ , M} 



(26f) 



where the numbers T2k (Tab. [D are defined by OlQct and satisfy the recurrence OlOct . 
Proof. Computing the integral in ( I26bl ) gives 

;,,(x;x,.,Ax) = J _ Z/'- ^"Z^^ ^ '2) -L^A^''^) 



m=0 
M 



ZCh^ 'V' + 1 1 /-^ - -^z Y 

, m + 1 ^\ n II Ax / 

m=0 \n=0 ^ ' 

M (mi , \ 

m=0 V/i=0 ^ ' 



m + \\(x-Xi\' 1 1 -(-1> 



m+l-n 



m + l\/x-x,V' 1 1 



(27) 
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where in the last Une of ( l27l l Yjn=o was changed to Tj'n=o because n - m+l - 
This proves that both ph(x; x,, Ax) and p /(x; x,, Ax) are of degree M. Since 



i-(-iy 



m+\-n 



m + \- n-2k+\ e No 
m+l - n - 2k A; e No 
0<n-m-2k<m ke No 

upon substituting 2k :- m- n, ( l27l ) becomes 



1 _(_iy"+i-" 
1 _(_iy"+i-" 

Q<2k< m 



m 

0<k<V-\ 



= 1- (-1)0 = 0. 

(28a) 
(28b) 
(28c) 



m=0 V't=0 ^ ' 

and, using ( IA.3I I. ( |29] > reads 



(^) 



= ZZ 

m=0 A:=0 



m + l 



22* (m + - 2;t 



y - \ m-2k 



I Ax / 



(29) 



M A^J 



/?/(x;x,-, Ax) = ^ ^ 



22* (^" + 2;!:+ 1) 



^ + 2A: + 1 



I Ax / 



f=0 ^^ i:=0 

Using the identity iAAi and changing the summation index £ to m gives 



(30) 



/7/(x;x;, Ax) = ^ ^ 



fll=0 /(:=0 



22*(2A:+1)\ 2k 



m + 2k' 



(31) 



which proves ( I26dl >. In practice, the coefficients cy;„ are computed by solving a Vandermonde sys- 
tem ||29], and the linear system ( I26dl ) must be solved to compute the coefficients c/,^ |0]. The general 
solution can be obtained using backward substitution without making reference to the basic recon- 
struction relations This alternative, matrix-algebra-oriented, proof of Lemma [3TT] is given in i l3.2[ 
However, the solution can be obtained immediately, by observing ( l26eD that the relation between 
Cf^ml and Ch,„^2k^m+2k)\ in ( I26dl ) is identical to the relation between f "'\x) and Ax^''h^'"'^^''\x) in OlOaH . with 
the only difference that the upper limit of the sum is finite. The inverse relation is exactly analogous 
to dlObl l. because, using ( l26eD in the right-hand-side of ( I26fl l 



— r y, T2s Cf^^,^ (m + 2s)l = — V t2s 



m. 

.9=0 
l^if k 

k=0 \s=0 



s=0 



/I it:!" I J \ 

(m + 2.9 + 2^)! c,,„,^,,^,, 



22f (2/"+ 1)! 



22*-2' (2;t-2i+ 1)! 



(m + 2k) ! c/,„,^2i. 



(m + 2A:) ! Ch^ 



s=0 



T2s 



where we used ( fl5l ). and ( IA.3I I and ( IA.2I ). This completes the proof 



(m + 2^ + 2^)! CA„„,„,, 
22^ (2^+ 1)! m! 

(32) 

□ 



The extension of the above results (Lemma 13.11 1 to infinite power-series (assuming that they are 
convergent) is straightforward. 

3.2. Matrix inversion proof of Lemma [3. II 

In this section we summarize the matrix inversion relations which can be used for an alternative, 
matrix-algebra-oriented, proof (Lemma 13.41 * of Lemma 13.11 By ( I26dl ) the coefficients cy;, of p / are ex- 
pressed as linear combinations of the coefficients c/,„ of pi,. This system ( I26dl >. whose solution expresses 
Ch,, as linear combinations of cy„ is the deconvolution linear system 12, (3.13b), p. 244]|f| Since the sum- 
mation relations ( I26dl > involve increments with step 2, we can split ( l26dD into 2 independent linear 



^more precisely, the system in 0, (3.13b), p. 244] relates n\cf^ with hi!c;i, 
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systems 



i=0 



for [c. 



^ChJ' ( I33al > and for [q, 
1 



Z '-/iM-l-2f+2t / 



2;t 

M - 1 - 2^ + 2/t' 



2fc 



.=0,...,Lfj 



(33a) 
(33b) 



( I33bl >. respectively. In matrix-form, we have 



1 

■•■ 
■■• 



1 lM-2 
(2+ 1)22\ 2 

1 





1 IM 

(4+ 1)241,4 

1 /M' 

(2+1)242 
1 





'^''m-21#J 




•^4-21 # J 












'-'lM-2 




C&-2 











(34a) 



(Lf J,M) 



1 

■■■ 
■■■ 



1 



M- 1 
(4+l)24\ 4 
1 lM-\ 

(2+ 1)24 2 
1 



CllM- 



'^/m-1-2iVJ 



(34b) 



t^(L*^J,M-l) 



where the matrices f/f^Mj^^) ( I34al l and ( I34bl l are upper unitriangular Il30ll . The corresponding 

linear systems ( l34l ) can be solved using backward-substitution lf30ll . To obtain the general solution, we 
initially remind, without going into the details of a formal proof, a standard result of matrix calcu- 
lus [30], concerning the inverse of an upper unitriangular matrix. 

Lemma 3.2 (Inverse of an upper unitriangular matrix). Let U e M"^" be an upper unitriangular 
matrix 



Uij -1 1 < i < n 
■' \ < I <n 



1 Ml,2 

1 



U = 








M2,«-I M2,n 



1 



(35a) 
(35b) 

(35c) 
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Its inverse U ' exists and is also an upper unitriangular matrix 



Uij —I I < i < n 
Uij -0 



J < ' 
I < i < n 



1 Ml,2 

1 







U2.n-l U2.n 



1 





^n-\,n 

1 



(36a) 
(36b) 

(36c) 



whose nonzero elements Uij (j > i) satisfy the recurrence relations 



— ^ ] l^n-k.n-k+t ^n-k+C,n- 



k+s 



I <k <n 
I < s <k 



(36d) 
(36e) 



Proof. It is straightforward to show, by induction, that detC/ - 1. The proof by induction of ( l36l l is 
a simple exercise of matrix calculus, directly obtained from the backward-substitution algorithm for 
solving Ux^b [30]. □ 

This recurrence is applied to compute the inverse of the upper unitriangular matrices 0341 1 of the 
linear system ( I26dl ) of Lemma 13 .11 



Lemma 3.3 (Inverse of the matrices in Lemma I3.1D . Assume N < LyJ + 1. Let U^n.m) ^ 

upper unitriangular matrix whose elements are given by 

Q<k<N -l-€ 



be an 



{U{N,M))N-t.N-t-k - 
{U{N,M))N-t.N-t - 1 

{U{N,M))N-t.N-t+k - 



M . 



1 



{2k + 1)2 



Ik 



M-2{ + 2k 
2k 



0<k< 



; 0<^<A^-1 ; A?<LyJ + l 



(37a) 



Its inverse t/, J is also an upper unitriangular matrix whose elements are given by 



{N,M) 

N 





1 



0<k<N-l-£ 



M 



(M-2{ + 2ky. 

(U(^f^^M)>N-t.N-t+k - T2k (^_2^)! — 0<k<£ 



; 0<{<N-1 ■ A?<LyJ + l (37b) 



where the numbers T2k (Tab. [T]l are defined by the recurrence ( llOcl l. 

Proof. To simplify notation let (U(N,M))ii - "0 ^iid (U^^j^-^)ij - uij By Lemma [3l2] ?/^^^^ is also an upper 
unitriangular matrix. It is easy to verify, by straightforward computation, using 

mi, that (3TU holds 

for < /" < 3. To prove that ( l37bD is valid for < f < - 1, by induction, suppose that ( I37bl ) is valid for 
l<(< m. Then, from ( l36e] > 



UN-{m+l),N-(m+\)+k 



■ ^ UN-{m+i)^M-(m+\)+s UN-(m+l)+s,N-{m+l)+k 
s=l 

k 

■ ^ UM-{m+\),N-(m+l)+s UM-(m+l-s),N-{m+l-s)+(k-s) 



(38a) 
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and since 5 > 1 m + I - s < m, we may replace UN-(m+i-s),N-(m+i-s)+{k-s) in ( l38aD by ( I37bl l, so that 

-1 /M - 2(m + 1) + 2s\ (M-2(m+ 1 -s) + 2(fe -.?))! 
«Ar-(™.i),A.-(,„.iH*-2,22. (2, + 1)1^ 2s r'-'' (M-2(m+l-.))! 



it 

/it ^ 



-T2,t-2. (M-2(m+l) + 2ky. 
22' (2.5+1)! (M-2(m+l))! 



Z 



-Tlk-ls 



-f 22^ (2i + 1)! 
1 y 



(M-2(m+ l) + 2fe)! (M-2(m+ 1) + 2A:)! 

= T2k (38b) 

(M-2(ot+1))! (M-2(m+l))! 



because, setting { :- k - s 



k-l 



(2,9+1)! f;^ 22*-2f {2k-2{+\)\ ^ ^ 

by OlOcD . This completes the proof of ( I37bl l by induction. □ 

Lemma 3.4 (Solution of the linear system ( I26dl l). The solution of the linear system ( I26dl l is given 
bv (l26B. 



Proof. The unitriangular matrices U^^_M^ l^^^ ( l34aD and f/(LM^j^_i) ( I34bl l are of the type defined in 

Lemma [3^ Using the result ( I37bl ) of Lemma 13.31 for the inverse matrices f/,7L, and J/tL , , , 
(Lf-J.Af) (L^^J.^w-l) 

the solution of the linear systems ( I34l > is 

(M-2^ + 2;t)! M 
Chu-ic = 2j '"2*'^/«-='+2* {M-2€)\ ^ =0. • • • ' Ly J (39a) 

■A (M- 1 -2^ + 2A:)! M-1 , 

where the numbers Tik (Tab. [D are defined by the recurrence ( llOcD . Since 

— 

m^M-2e ^21 =M - m^^=L — ^ — J (40a) 

M — m 

M -2e -l^2( +\^M -m^(^[_ — - — J (40b) 
the 2 solutions ( [39l l can be grouped into ( I26f] l. which completes the proof □ 

4. Error of polynomial reconstruction 

We consider in this paper reconstruction on a homogeneous grid (recall that ^ hold iff Ax = const). 
The reconstruction polynomials are computed by interpolating f{x) sampled on an appropriately cho- 
sen stencil (Definition l4.1D . We examine the relations and order-of-accuracy of polynomial reconstruc- 
tion (Definition l2.3[l on an arbitrary stencil S,,m_,m+ (Definition l4.1D defined on a homogeneous grid. The 



WENO i5l.^. lTlIll3l . l3 . [l5ll schemes are based on the convex combination of polynomial reconstructions 
on a family of substencils. For the development of the order-of-accuracy relations, it is necessary to 
develop results on the approximation-error of polynomial reconstruction for the general stencil S,,m_,m+ , 
around point / (not necessarily contained in the stencil), with M_ neighbours on the left, and M+ neigh- 
bours on the right (Definition l4.1D . 
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4.1. Polynomial reconstruction 

The part concerning the approximation of fix) by a polynomial pf{x; S,;m_,m+, Ax) is found in most 
textbooks of numerical analysis ll20L[2lll . It is only briefly included here for use in deriving the results 
concerning the approximation of h{x) by the polynomial ph(x\ S;,m_,m+, Ax) which forms a reconstruction 
pair with p / (Definition l2.1l >. To obtain the relations concerning p7,(x; S,;m_,m+, Ax) it is not very practical 
to work with the Newton divided-differences form of p/ ll20l . I21J I. which are widely used in WENO 
theory 0,I2,SI3, 4]. It is, instead, preferable to work with the standard form of/?/ expanded in powers 
of (x - X,), whose coefficients can be readily expressed (Proposition I4.5D from the coefficients of the 
inverse of the Vandermonde matrix jsills^ corresponding to the stencil S,,m_,m+ (Definition 14. 11 1. This 
representation of p / allows direct use of the formulas relating the coefficients of ph and p / (Lemma lS.lD . 



Definition 4.1 (Stencil). Consider a 1-D homogeneous computational mesh 

Xi - xi+ (i - l)Ax Ax = const e M>o (41a) 

Assume 

M M_ + M+ > (41b) 

The set of contiguous points 



5/,M_,M+ •- 



M_,--- ,! + M+} (41c) 



is defined as the discretization-stencil in the neighbourhood of /, with M_ neighbours to the left and 
M+ neighbours to the right. The stencil S,;M_,Mt Hid contains M + 1 > points and has a length of M 
intervals. If M± > then the stencil contains the pivot-point /. If < then the stencil does not 

contain the pivot-point /. We will note 

[Si,M_.Mj [Xi-M_,Xi+Mj c M (41d) 

the interval defined by the extreme points of the stencil. □ 

Remark 4.2 (Stencils and notation). In our notation the stencil is defined by a reference (pivot) 
point and by the number of neighbours M± on each side of point / (Definition 14.11 1. The position of 
the pivot point / in the stencil is arbitrary. This is necessary for obtaining relations for all of the WENO 
stencils with reference to the same point /. In the following developments, there appear quantities 
depending both on M± and on / (and eventually on the values of / sampled at the points of the stencil). 
We will systematically note these quantities as functions of the stencil S,,m_,m^- On the other hand, 
there appear quantities, which depend on M± but not on the pivot point / (neither on the values of / 
sampled at the points of the stencil). We will systematically note these quantities as functions of M_ 
and M+, and not of S,;m_,m+- This difference is important when considering order-of-accuracy relations 
(eg Corollary [4]9l). □ 

Definition 4.3 (Vandermonde matrix on S, m_ m+ )• Let M : = M_ + M+ and assume M > 0. The matrix 

;;^-y e R(«+i)x(«+') with elements (^^-y),; 



{-MJf (-M_)i ■■■ (-M_)*' 

M := M_ + M+ > (42) 

(+M+)° (+M+)i ■■■ (+M+)'^ 



is the Vandermonde matrix llsil [3^ defined on the stencil S,,m_,m+ (Definition iHl. Since is a 
Vandermonde matrix, its inverse '^*V~^ exists Hi, HI. The elements of ^"V ' e rCw+Dxcm+d 
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Lemma 4.4 (Inverse Vandermonde matrix on S,;m_,m+)- Assume the conditions of Definition 14.31 
Then the entries of the inverse of the Vandermonde matrix '^* V ( [42] l on S,,m_,m+ are given by 



M+l- 



n=0 



n + i - I 
n 



(qV )i+n,i 



V/,je {I,-- - ,M+ 1} 



(43a) 



where is the inverse of the Vandermonde matrix on S^\m = {/,■■■,;+ M] (Definition l4.3l l. whose 
entries are given h"^ 



(2^y-')„ = (-iy-^£ 





A: - 1 




/-I 



V/, j e {I,-- - ,M+ 1} 



(43b) 



Define 



1^ 



(43c) 



T/ien ^/le following identities hold 



'm+l,C+M. 



< ;t < M 
< m < M 



(43d) 



M 



6 {-M_, ■ 



(43e) 



Proof. Since ^* V ( |42] | is an (M + 1) x (M + 1) Vandermonde matrix on M + 1 distinct nodes its inverse 

J^+y ' exists ll29l[33ll . Macon and Spitzbart ll33.l29l1 have given explicit expressions for the inverse of the 
Vandermonde matrix on integer nodes. To prove ( I43bb we start from llssl. Theorem 1, p. 973], giving the 
inverse of the Vandermonde matrix on n equidistant nodes on [0, 1], ie on (n - l)xi = (/- 1) V / e {1, ■ • ■ , «), 
as 



1 



n-1 





k-l 




/-I 



(44a) 



are the unsigned Stirling numbers of the first kind l22Ll23Ll33ll satisfying 



n + 1 


=n 


n 


+ 


n 


k 




k 




k - 1 



n - k 
m + i ~ k 



m- If 



« - 1 
/t- 1 



=(-!)" («-!)! 



n~ k 

m - 1 
n - 1 



* Proof of j43d) is most easily obtained using Proposition l4.5l and proof of <43e) is most easily obtained using Proposition l4.6l 
which are proved below. Notice that j43d) is not used in the proof of Proposition l4.5l nor is <43e) in the proof of Proposition l4.6l 
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which directly impHes, setting n - M + I, 



i-1 



, /,;6{1,-- - ,M+11 



M+I 



= (-1)-' J] — - 



k=\ 





k- 1 




i-1 



(44b) 



Obviously, M' ' and M' ' in ( I44bl ) are scaling factors (for M + 1 equidistant nodes on [0, 1] we have 
M Ax = 1). This is clearly seen by writing the Vandermonde matrix on S,;o,m < I42D as 



^V:^[(i-iy-\ ijeih--- ,M+1]] 



i-1 



e {I,-- - ,M + 1) 



[m^"' dej, {J e{l,--- ,M+ 1)] 



and since [m^ ' 5f y e { 1 , ■ • ■ , M + 1 }j is a diagonal matrix 



ijell. 



,M+1] 



€-1 
M 



i-1 



(,je{l,--- ,M+ 11 



(44c) 



(44d) 



which, by ( I44bl l. proves ( I43bl ). 

To obtain the final expression ( l43aD . we observe that, for M := M_ + M+, the stencils S,,m_,m+ (corre- 
sponding Vandermonde matrix jJJ^ V; Definition [4j3) and S,-m_,o,m (corresponding Vandermonde matrix 

; Definition l4.3l l correspond by Definition l4.1l to the same set of points {/ - M_, • • • , ; + M+). Therefore, 
], by the uniqueness of the Lagrange interpolating polynomial lH^], we have (using 
the notation of Proposition [4l5] l 



Vxe: 



P/(^; S,;m_,m+, Ax) - /?/(x; S,_M_,o,M, Ax) \i f ^ r\ ^ 

" / fc LXi-M_ 5 Xi+M^ 1 



(44e) 



the only difference being in the choice of the pivot point (x, for S,,m_,m+ and x,_m_ = x, -M_ Ax for S,_m_,o,m) 
used for the representation ( I45bl l of the interpolating polynomial of /(x) on the nodes {/-M_, ■ ■ • , ; + M+}. 
By ( I45bl >. ( I44el l reads 



M 



M 



M 



Z c,,s_,„,„, + M_)- = ^ X c,,s,_„_,„.„,.,p(M_)« (^) 

^■=0 5=0 77=0 ^ ' 

M rM-m / \ \ m 



implying 



m=0 V n=0 



M~m 



VxeM 



n=0 



Im + n 



i(M_)" V OT e {0, ■ ■ ■ , M} 



(44f) 



(44g) 



which by ( I49al > gives, V / e C[xi^M^,Xi+Mj 



Z (J^*^ )m+I,f+M_ + I fi+e - Z ^(o'^ )Mj+n+l,.s+l fi-M_+s 



«=0 Vj=0 
M M-m 



)m+n+\,s+\ fi-M_+s 





m + n\ 




[ n ] 



s=Q n=0 
f=-M_ V ,!=() \ I 



m + n 
n 



(M_)" 
{MJf 



m+«+l,f+M_ + l 



fi+( 



(44h) 
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and since fi+i {t e {-M_, • ■ ■ , M+}) are linearly independent we have 



\M- )m+\X+M- + \ 



+ „My-K Vm€{0,---,M) 

tt=0 ^ ' 

which proves ( I43al >. 

To prove the identities containing VM_,M^,m.k ( l43cD . notice that the elements of ( l42] l read 

fi:v)„=„-i-«_y- ,44,) 

Explicit expression of the elements of the product (|JJ" V ') ■ (JJ* V) = /m+i (where /m+i e M('^+i)x(*'+i) is the 
identity matrix) yields 

- [cr') ■ (:::v)L,,^. = v._,.,,„„ (44k) 

and as a consequence ( I43dl ). To prove ( l43eD . consider the error dSlbl l of the polynomial interpolation 
Pf{xi + ^Ax; S,_M_,M+, Ax) on the stencil S,-,m_,m+ (Proposition l4.6l l. By construction, we have 

Pf{xi + Mx; S,-,M_,M+, Ax) = ^ ^/(x,- + thx; Si,M-^M,,Ax) = e {-M_, • ■ • , M+] (441) 

which, using dSlel l and ( |51g[ ) in dSlbl l. proves ( l43eD . □ 
Proposition 4.5 (Lagrange polynomial reconstruction on S, m ,m+). Let 

M 

\~\ j X — X \'" 

Ph{x\ Si^M__M^,Ax) :=2^ c/,.s,„_,„^,m ( ) (45a) 

m=0 

v"' /X — x- V" 

/7/(x; S;,M_,M„ Ax) :=2_, c/,s,„_,„,,m (^^^ j (45b) 

be 2 polynomials of degree 

M:^M^+M+ (45c) 

constituting a polynomial (Lemma l3.1l > reconstruction pair (Definition l2.1D ph - R(\-Ax){pf)- Assume that 
the polynomial Pf(x; S,;m-,m+, Ax) is obtained by interpolation of the values of /(x) on the points of the 
stencil S,_m_,m+ (Definition l4.1l l. Then 

/7;,(x,- + ^Ax;S,;M_,Mt,Ax) = ^ ah,M-,M^,M)fi+e (45d) 

f=-M- 

/7/(x,- + ^Ax;S,;m_,m^,Ax) = ^ af^M_^M^,e(^)fi+e (45e) 

where ah,M-,M*,M) Of^M-,M+,M) '^^^ polynomials of degree M in 

^ := ^ (45f) 

Ax 

with coefficients depending only on the 3 indices (M^, M+, C) 



r (45g) 



V V '^2k{m + 2ky -I 

CthM-,M+A^) ■- 2_i 2-i i ^f^- )m+2k-¥\,e+M_+\ 

M 

af,M^,M.AO := Yj^M:v'')m+i.e+M^^i f (45h) 

m=0 
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where Qv' ) ij are the elements of the i nverse Vandermonde matri x on S,_m_,m+ (Lemma and the 
numbers T2k (Tab. [Hi are defined by dlQcl l and satisfy the recurrence OlOcl l. 

Proof. Define 

Xi+c + tt^x - < { < M+ (46a) 
/(^,+f) -M_<^<M+ (46b) 

The M+ 1 coefficients c /,(s/,„_,m+ ),m («^ = 0, ■ ■ ■ , M) are computed by equating the polynomial p /(x,+f ; S/,m_,m+ , A^) HBhi 
to known values f^e 

(47) 

fi+M^ - Pf(Xi-M+', Si^M_,M^, ^X) 

Expanding ( I47l > results in an (M + 1) x (M + 1) Vandermonde (Definition l4.3l > linear system 



(-M_)° (-M_)i ■■• (-M_)*' 
(+M+)° (+M+)i ■■■ (+M+)'^ 



^f,Si.M-.M+fi 
'^f,Si,M-.M+'M 



fi-M_ 

fi+M+ 



(48) 



y 



Hence (Definition |4]3) 

'^/,S,,M_.M+ ,m - ^ (m!^ )m+l/+M_ + l /i' 



Vm e {0, ■ • ■ , M) 



1 - ^ 



Ca,s,„ „^,m =— : y, T2* c/,s,„ „ ,,„+2/t (m + 2ky. Vm 6 {0, ■ • ■ ,M) 



(49a) 



(49b) 



k=0 



where we used the deconvolution formula ( I26f] > for computing c/,_s,„_„^,m- Injecting ( l49aD into Hdhi we 
have 



M+ / M 



„ V 'y '^ ^M, ,.-1, 

m=0 Vf=-M_ y fc-M_ Vmj=0 

proving 045eD and ( I45hl l. Injecting ( I49bl l into ( l45aD we have 



./i+f (50a) 



M 



/?/,(.x; + ^hx; Si,M-,M* , Ajc) = ^ ^ 

m=Q \ k=0 

M L— J— J , ^ , M 

^ T2A:(m + 2fe)! 



T2t(w + 2ky. 



m=Ol yt=0 



Z (m!^ )m+2*:+l,f+M_+l/i+ 



v-l v-l v-l T2t(m + 2fc)! ^ -1 
Zj Zj Zj ^M^*^ )m+2*:+I,f+M_ + I 



f=-M_ \ m=0 \ k=0 



i+C 



proving ( I45dl l and (|45g 



(50b) 

□ 
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4.2. Approximation error of Lagrange polynomial reconstruction 

Of course the accuracy relations for the approximation of f{x) are well-known ll20ll . but this section 
( ii4.2D is concerned with the accuracy of the approximation of h{x), using Lagrange polynomial recon- 
struction based on the knowledge of the values of f(x) on an arbitrary stencil defined on a homogeneous 
grid (SUl. 

Proposition 4.6 (Error of Lagrange polynomial reconstruction on S,m Let Pf(x; S,;m ,m+, Ax) 

and pi,(x;Si^M-,M^,^x) be a polynomial (Lemma I3.1D reconstruction pair (Definition 12.11 1 p/, - R(\-Ax)(Pf), 
satisfying the conditions of Proposition |4j5] Then, pjix; S,;m_,m+, Ax) approximates f{x) to 0(Ax*^^'), and 
Ph(x; S,,M_,M+, Ax) approximates h(x) to 0(Ax'^'^') 



Ph(x; S;,M_,M+,Ax) ^h(x)+Eh(x; S;,m_,m,,Ax) ^h(x)+0(Ax^^^) 
Pf(x; S,;m_,m+,Ax) =/(x)+£'/(x; S;,m_,m,,Ax) =/(x)+C>(Ax'^+') 



(51a) 
(51b) 



where the approximation errors constitute a reconstruction pair E/, = R(\;Ax)(Ef) (Definition 12.11 1 and, 
VA^Tj > M + I, are given by {assuming f and h are of class C'^'tj+I) 



Ehixi + ^Ax; S,-,M_,M,, Ax) = J] ///,,m_.m,..(^)Ax'/;.^"^ + 0(Ax^"+>) 



s=M+\ 



Ax'hf^ + 0(Ax™--^') 



(2{+ 1)! 

Efixi + ^Ax; S,;m_,m,, Ax) = J] ///,m_,m„.(^)Ax^/;^'' + 0(Ax^™^') 



s=M+\ 



where i-ih,M_,M+.M) o^c? ^J■f,M_.M^,.M) are polynomials of degree s in ^ (45 



(51c) 
(51d) 
(51e) 



*:=0 
1 



{s-2k)V 

M 



T2kVM_,M+,m+2k,s- 



{m + 2k)\ 



m=ay k=o 



\ m=0 



(51f) 



(51g) 



where VM_,M^,m.s are defined by ( I43cl >. and the numbers T2k (Tab. [T) are defined by ( ll9cD and satisfy the 
recurrence OlOcH . 

Proof. To prove idlhi we start by Taylor-expanding /j+f in ( l49aD . and using H3di 



/'A', 



''/■S,-.M_.M+,m — ^ (m*^ 



m+l/+M_+l 



A't., /' M+ 



V J=0 



^ Ax'/ '' 



.5=0 Vfc-M_ 



+ C>(Ax^"+') 



Ax'' f*^'* ''^ Ax'' f*'''* ''"'^ Ax* f '■■^ 

2 VM_,M„,„,. + 0(Ax^™+i) = 2 ^ + 2 ,„ + C»(Ax^-+') 



.5=0 



j=0 



s=M+l 



Ax"'/. 



<m) Nt.j 



Ax'/.' 



is) 



+ C»(Ax^"+') 



(52) 
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Injecting ( [52l l into ( I45bl l. and replacing /(x, + ^Ajc) by its Taylor-polynomial, we have 

Ef(xi + ^Ax; S;,M_,A/^., Ax) =p /(x,- + ^Ax; S/,m_,m+, Ax) - /(x,- + ^Ax) 



m=0 
M 



Ax'"/;.' 



.(mj) W. 



Ax* f 

VM_,M,.m,.— + 0(AX^" + 1) 



.s=M+l 



Vhi=0 

^ -Ax'"/' 



^r-/(x,+^Ax) 

■(m) M N. 



' m=0 s=M+ 1 



r-/(^,+^Ax) 

Ax' f*'' 

VM_,M„«,,.. + 0(Ax^"+i) 



\ M Ntj a s f (s) 

^ ^ .. /-m , ^/a..A't., + 1 



A'tj a ..j As) Ntj ( M 



Ax^f'*' 

r+Yj Z VM_,M„m,.— + 0(AX^™^') 



^ £ -Ax*./;: 



.!=M+1 



.s=M+l 



Vm=0 



Ax*/.' 



s! 



+ C>(Ax^"+') 



proving dSlbl l. OSleD and (|51g 



To prove dSlal l we use the expression ( l52] l for c/,s,„ „ ,m in ( I49bl > to obtain 



c/i,s,« „^,m =— : y. {m + 2k)\ 



k=0 



{m + 2ky. 



Ax*/.' 



.m+lk Am+2k) L^J 



Ax'"+2*/ 



yt=0 



ot! 



i=0 



T2A:(m + 2k) 



s=M+\ 



• C>(Ax^™+') 



2 

V.s=M+l 



■''M_,M+,m+2j(:,.s- 



Ax*/*^^ 



+ C>(Ax'^"+') (54) 



Injecting ( [54l ) into ( I45al l. and replacing /z(x, + ^Ax) by its Taylor-polynomial ( fl6l >. we hav^ 



jr(m+2k) 



E/^x; + ^Ax; S/,M_,M+, Ax) = /7/,(x,- + ^Ax; S,-,m_,m+, Ax) - /i(x,- + ^Ax) = /-i Zj ^^'^ ^J^f ^ 



m=0 Jl=0 



-ZEE 



A'tj J 



/ , i;m Av-'f^*-* «Tj J f. m+2k Am+2k) 

T2k(m + 2ky. Ax/. ^ T2t Ax /. 

; yM-,M+.m+2k.s j g - / , / , j 

m=0 k=0 s=M+l ' ' m=0 k=0 



i"" + 0(Ax'^™+') (55a) 



which simplifies to 



Ax) = 2 



m=0 k=Q .s=M+I 



T2/t(m + 2^)! 

; '^M- ,M+ ,m+2k,s T 

ml si 



Ax*/: 



s A-') 



M L^J 



Zj Zj «,,! 



m=M+I *:={) 



m=0 t=Li£^J+i 

Using ( IA.3I I and ( IA.2I I. ( I55bl > reads (the summation indices on line 1 remaining unchanged) 



(53) 



(55b) 



£'/,(x,+^Ax;S;,M_,M+,Ax) = Z Z Z 



/v„ (' LfJ 

Z 



i=M+l 



-T"2* _^.?-2y(: 



« («-2fe)! 



.s=M+l \m=0 k=U 



T2k(m + 2k)] 
si ml 



VM-,M+,m+2k,s 



Ax*/^-'^ 



Ax*/.w+ y y ^ii^^s-2k 



Ax*/*' + 0(Ax^"+i) (55c) 



JVtj 



LfJ 



(s-iey. 



Ar' /-"W + 0(A.r"-+') '"^-'' f V rzr A.r"-^^^ /<"'+^^'(.v) _^ o(Ax"™-^') 
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and defining Hh,M-.M+.M) by dSlfl l we obtain dSlaH and dSlcH . 
Finally, using dlOaD in dSlcD 



f^h,M-,M„s-.-,^ Z-i (-2/+ IV 



.5=M+1 ^=0 



(2{+ 1)! 



which, by dA.3l l and dA.2l ). proves dSldl l. 



(55d) 



□ 



Proposition 4.7 (Approximation error of Lagrange polynomial reconstruction on S,,m ,m+)- Asssume 
the conditions and definitions of Proposition 14. 61 Then 



Ax) = J] i/,,M_,M„«(^) Ax" h^^Kxi + f Ax) + 0(Ax^"+') 

n=M+\ 

A^TJ 

£/(xi + ^Ax; S,,M_,M„ Ax) = ^ ^/,M_,M„„(a Ax" f"\xi + ^Ax) + 0(Ax^"+i) 



(56a) 
(56b) 



n=M+\ 

where /l/,,M_,M+,n(^) and Af,M-,M+.n{0 are polynomials of degree n in ^ (45 

n-M-l , 



f=0 
n-M-\ 



(~^) 



(56c) 
(56d) 



where fJ-h,M-,M+,iM) '■^ defined by d51fl l and jj.f^M-,M^,n(^) *s defined by ( |51g 



Proof. Taylor-expanding /; in d51eD . around the point x, + ^Axjf] we have 

J i3LA//(-^)(x,.+^Ax) 



^/(X; + ^Ax; S;,M_,M+, Ax) = ^ yU/.M_,M+,s(^)Ax' 

s=M+l 



C»(Ax^"+') 



2 Z /■"''(-^' + ^^^) + 0(Ax'^"-') 



.s=M+l f=0 
Ntj II- M- I 



= Z Z /^/.M-,M„«-f(^)^ Ax« /("'(X,- + ^Ax) + 0(Ax^"+') (57a) 



n=M+l e=0 



which proves d56bl l. 



Ajc' /<"+'"' (jc + f A,v) + 0(Aa-""+1) 
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Replacing /j.**^ in dSlcH by its expansiorQ in terms of the derivatives Ajc^ + ^Ax) we have 

Ehixi + ^Ax; S,,M_,M, , Ax) = V ;U,,,m_,m„.©Ax' V J— i— ((f - 1)^+' - + 1)^+') A/ h^-'^'\xi + ^Ax) 



.s=M+I 

+C>(Ax'^"+') 
+0(Ax^"+') 

A'tj n-M-l , , , 



= Z Z /^'>.M-.M,,.(^)|^^ ((^ - i)'"' - + i)'"') Ax« + ^Ax) 



n=M+l e=0 

+C»(Ax^"+') 



(^ + 1) 



(57b) 



which proves dSGal l. Obviously, by ( I56dl l. deg(/i/_M_,M+,n(f)) = «■ It is easjo to verify that, by ( l56cD . 
deg(/i/,,M_,M+,n(^)) = n, which completes the proof □ 

4.3. Approximation error ofh^^t, and of fl 

One of the principal uses of the reconstructing polynomial being the numerical approximation 
of /j' :- f'(xi) via we give in this section the relations concerning the approximation error of 
^i±i, '■- ^(^i ± 5^-*^) (Corollarv l4.8l > and of /j' (Corollarv l4.9D . which are readily obtained by application of 
Proposition l4.7l 

Corollary 4.8 (Accuracy at / + 5 of Lagrange polynomial reconstruction on S, m ,m+ )• Let p fix; Si^m-.m^ , Ax) 

and p/^x; S,,M_,M+, Ax) be a polynomial (Lemma I3.1D reconstruction pair (Definition l2.lt p/, = R{\-a.x)(p f), 
satisfying the conditions of Proposition l4.5l Then, the reconstructed value at Xj^\_ :- x,- + ^Ax, which will 

be noted h^^^ ^ approximates :- h{x^^i_) to 0(Ax*'+') with M :- + M+ > 0. The error of the ap- 
proximation can be expanded in powers of Ax with coefficients involving the derivatives h^"'\ :- /!*'"^(x;^ 1 ) 



'*s,-,„_.„^,i+i :-/?/,(-X/+^; S,-M_,M+, Ax) 



- Z '^M_,M^,e fi+e 

t=-M- 

A'tj 

1 + Yj Am_,m„.Ax^/z|;', + 0(A./'™+i) 



(58a) 
(58b) 

(58c) 



s=M+\ 



0(Ax^+i) 



^ Approximating h{() (which was assumed to be of class C*^ in Lemma I2.2> in {6a) by the corresponding Taylor-polynomial 
(Taylor-jet) of order A'tj 1 27, pp. 219-232] around ^ = .v + ^hx yields, VA'tj e N : A'tj < A', 



f(x)=l- p-^' (fy (L^_^/,m(x + fA.r)l + 0((f-.t-fA.vf"+') 



■•x+iAx (N-a 



-fAx)L^ + 0(Ax"--')=-l£ r - ^ 



dr] h'^^\x + f Ax) + 0(Aa*J +' ) 



"TJ \ 

Z ((^ - - + ^)'") A-^-' f''"^' + ^Ax) + 0(A.r~"+') 



k p [2) 



-E 



'((4r-(^)>E(f;>-i(4r-(^^r); 
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where the constants Am_,m^,s are given by 

s-M-l ^_^y 

Am_,M+,s '■= ^h,M_,M+,s(j) - 77 Tr^l^h,M_,M+,s-e(j) (58d) 

(t + l). 



with Ah,M_.M^,.M) being the degree s in ^ polynomial defined by ( l56cD . Hh,M-,M+,si.O being the degree s in ^ 
polynomial defined by dSlfl l, anC^ 

om^.m^/ ■- o:h,M-,M+A\) (58e) 
with ah,M-.M+A^) being the degree M in ^ polynomial defined by ( |45g l. 



Proof. Using ( I45dl > and 056al >. in the definition oih^^^ ( l58aD . we have immediately 



f=-M_ 

Nil 

Yj ^,M_,M„.(i)Ax''M;\ +0(Ax^"+') (59) 



.s=M+I 



Eh{xi + \ Ax; S,-m_,m^, Ax) 



and using the definition ( l56cD to compute /1/,,m_,m+,.s(5) completes the proof. □ 

Corollary 4.9 (Order-of-accuracy of Lagrange polynomial reconstruction). Assume the condi- 
tions of Proposition |4j5] Then 

h ■ 1 -h 1 '^Tj 

. '^'-■"-"-'-^ =/;■'+ Z AM_,M„„Ax"/i("^'^ + 0(Ax^--') = + 0(Ax«->) (60) 

where /isi_,„ - '^s, ,„ „ OSSal l. and ^/le constants A.M_,M^.n are defined by dSSdl l. 

Proof. The constants Am_,m+,« dSSdl ) depend only on the 3 indices iM-,M+,n), and not on the point 
index / (Remark 14.21 1. because the polynomials i^h,M-,M+.s-e(^) dSlfl l are also independent of the point 
index Hence, we have, by d58ct , 

k i I : ,^1 1 + y Am M,M'k'\ +0(Ax^"+') (61) 



s=M+l 



Subtracting d6T] ) from dSScD yields 



+ 2j ^M_,M„..Ar' — + OiAx^'"^^) (62) 



Ax Ax ^—^ Ax 

s=M+ 1 



and using the exact relations (HJl we obtain d60l ). □ 



^ Notice that Shu |3], following a different route, has shown that 

M+l M+1 

2 W (M_-?+l) 

p i^m q ^ m 
1^' 5 * p 

aM-,M,.f - 2j 

m=e+M-+\ , . 

j j (m-p) 

p = 
p i^m 

is an equivalent expression for the coefficients au-M^f iHD- 
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Remark 4.10 (Order-of-accuracy). The previous result (Corollarv l4.9l l illustrates that the 0{l^x'^*^) 

accuracy in approximating /' is achieved, using Oil^x'^* ' ) interpolates for /, because of the exact recon- 
struction relations (Lemma 12.21 1 . Liu et al. note this as an (9(Ax'^) accuracy increased to (9(Ax'^+') 
at one chosen point, viz jc,. □ 



5. Interpolating and reconstructing polynomial 

We briefly summarize how the existence and uniqueness properties of the interpolating polynomial 
carry on to the reconstructing polynomial. Consider first the general case of a polynomial reconstruc- 
tion pair ( ij3.lt . Combining the existence (Lemma l3.lt and uniqueness (Remark 12. 8D of polynomial 
reconstruction pairs, we can formulate 

Theorem 5.1 (Vector spaces of polynomial reconstruction pairs). Consider the {M+ 1) -dimensional 
vector space of polynomials with real coefficients of degree < M in x, Mm[.jc]. Then the reconstruction 
mapping R(\-^x) (Definition l2.1l l is a bijection ofRmlx] onto itself 

Proof. By construction (Lemma 13.11 1 V p(x) e Rm[x] 3 q(x) - [R(1;Ax)(p)](x) e Rm[x], and inversely 
V q(x) e Rm[x] 3 p(x) - [Rj^,^^^{q)](x) e RmIx]. Furthermore, since the elements of Rm[-«] are continu- 
ous functions, the reconstruction pair q(x) - [R(i;&x)(p)](x) is unique (Remark l2.8l >. which completes the 
proof □ 

In his recent review of WENO schemes, Shu iQl stresses the difference between WEND interpola- 
tion and WENO reconstruction. In this sense, pf(x; S,-,m-,m+, Ax) in Proposition 14.51 is the interpolating 
polynomial of f(x) on S,-,m_,m+, and ph(x; S,;m_,m+, Ax) is the reconstructing polynomial (Definition l2.3D . Of 
course 

Proposition 5.2 (Lagrange reconstructing polynomial). Assume the conditions of Proposition l4.6l 

The Lagrange reconstructing polynomial ph(x; S,;m_,m+, Ax) approximates h(x) to (9(Ax*^^') but, unless fix) 
is a polynomial of degree < M, it does not interpolate h(x) on S,-,m_,m+, ie, if fix) is not a polynomial of 
degree < M, we have in general 

Phixi + Mx; S,-,M_,M+, Ax) 7i hixi + {Ax) 'it e {-M_, • ■ ■ , M+} (63) 

Proof. Proof is obtained by contradiction. It suffices to give an example where the inequalities 0631 * 
hold. Consider the reconstruction pair (Theorem l2.9l > 

/(x):=e^-^' ; ^x) = [/;(1;a,)(/)](x) = g,(Ax)e^-^' (64a) 

with gr defined by ( I19bl >. Consider the polynomial reconstruction of fix) (Proposition 14.51 * on S,.i,i. By 
( I45dl > and ( |45g[ ) 

PHix, + ^Ax; S,,i,i, Ax) = [\f - 1^ - ^) + _ f)^ + f._^ (1^2 + 1^ _ ^) (64b) 

We have f - 1 and f±i - 6=^^^', and evaluating p/,(x, + £Ax; S,,i,i, Ax) - /i(x,- + ^Ax), using ( I64bl l and ( I64al l. 
for ^ = -1, 0, 1, and for different values of Ax ieg Ax - ^), we verify (63^ . □ 

Most of the results of existence and uniqueness properties of the interpolating polynomial hold, 
with appropriate adjustments, for the reconstructing polynomial, because of Theorem 15. 11 We briefly 
summarize in the following those necessary to prove WENO reconstruction relations Ii3i,i^. 

Theorem 5.3 (Existence and uniqueness of the Lagrange reconstructing polynomial). Assume 
the conditions of Proposition l4.6l There exists a unique Lagrange reconstructing polynomial phix; S,;m_,m+, Ax) 
of the form ( I45dl > which approximates hix) to (9(Ax*^+'). 

Proof. Existence, with ah.M-.M^.M) given by ( |45g[ ), is proved in Proposition 14.51 bv construction. We 
know from approximation theory '21] that there is a unique Lagrange interpolating polynomial 
Pfix; S,;m_,m+, Ax) on S,-,m_,m+, and that the reconstruction pair phix; S,;m_,m+, Ax) = [^(1;Ax)(p/)](-^; S,'m_,m+, Ax) 
is unique (Remark l2.8D . which completes the proof □ 
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6. Examples of applications 



The analytical relations developed in the present work can prove quite useful in the analysis of 
practical WENO schemes, and more generally in the development of discretization schemes. Providing 
detailed analysis of such applications is beyond the scope of the present paper. We sketch, nonetheless, 
in the following, 3 applications (the complete proofs will be given elsewhere), to illustrate the useful- 
ness of the reconstruction pair concept, of the associated application of the deconvolution Lemma [231 
and of the explicit expressions for the Lagrange reconstructing polynomial. 

6.1. Representation of the Lagrange reconstructing polynomial by combination of substencils 

All WENO [3, 4] schemes for reconstruction on the general homogeneous stencil s,;m-,m+ (Defini- 
tion |4!T]| are based on the weighted combination of the reconstructions on + \ < M :- M_ + M+ 
substencilJ^ 

M± e Z : M = M_ + M+ > 2 
Qi,M--K,M,-KM. ^{i-M^+k,,--- ,M^-K, + k,} l<K,<M-l (65) 

with appropriate weights, which are nonlinear in the cell-averages f(x), to ensure monotonicity at 
discontinuities [35], and such that the weighted combination of the Lagrange reconstructing polyno- 
mials on the substencils, at regions where h{x) is smooth, approximates to 0(Ax'^'^^) [6] or higher l isll 
the Lagrange reconstructing polynomial on the big stencil S,;m _m , which (Proposition l4.6D is 0(Ax'^'^^)- 
accurate. The starting point for scheme design is the determination of the underlying linear scheme, 
ie the determination of weight-functions o-h^M_,M+,K„kM) ^ (0, ■■ ,^^s)), independent of fix), which 
combine the Lagrange reconstructing polynomials on the substencils exactly into the Lagrange recon- 
structing polynomial on the big stencil S,,m_,m+ 

Ph(xi + f Ax; S,;M_,Mt, Ax) = ^ crh.M_,M^,K„k,{^) Ph(xi + f Ax; S,;M_-t„M+-A:,+*:,, Ax) (66a) 

A-,=0 

Obviously, using ( I56al l in ( I66al l. the weight functions o-h.M-.M+,K„kS.^) must satisfy the consistency condi- 
tion 

^ 0-h,M_,M^,K„kM) = 1 (66b) 
k,=i) 

Shu [15] indicated examples of instances where it was impossible to find such weights, as well as in- 
stances where convexity of the combination was lost (presence of negative weights). The corresponding 
problem for the Lagrange interpolating polynomial 

Ax) (67a) 



is directly related to the Neville-Aitken algorithm ll2l]. pp. 11-13], which constructs the interpolat- 
ing polynomial on {/ - M_, ■ ■ ■ , / + M+] by recursive combination of the interpolating polynomials on 



Notice that the family of subdivisions <65> includes (when varying Kg) all possible subdivisions to substencils of equal length 
(M - A's intervals) whose union is the entire stencil S,_m_ . 
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substencils, with weight-functions which are also polynomials of x 11211. pp. 11-13]. Carlini et al. (Sf 
have given the explicit representation of the polynomial weights for the construction of the Lagrange 
interpolating polynomial on S,,m_,m+, by combination of the Lagrange interpolating polynomials on the 
K^ + l substencils, for a certain family of stencils/subdivisions. Liu et al. lll9n have extended the family 
of stencils/subdivisions studied. Since for the Lagrange interpolating polynomial case the weight- 
functions are polynomials, ( I67D are valid e K. 

The corresponding problem for the Lagrange reconstructing polynomial ( I66I I was studied, only very 
recently, by Liu et al. 1119(1 . It turns out that, for the reconstruction case, the weight-functions which 
satisfy ( l66l l are rational functions of ^, implying that dMl ) is valid € M \ So-,,j^_j^^j.^ , ie everywhere 
except at the union Sa-i,n,_ of the poles of the K^^ + 1 rational functions cr h,M-,M^,K,,kX^^ ^ fO' ■ ' ■ > ^s})- 
Liu et al. Ill9] studied the family of stencils S, |^mj j^_|^mj, for the K^^ - |^yj-level subdivision, in the range 
M e {2 ■ ■ ■ , 11}, and used symbolic computation to give explicit expressions of the weight-functions, and 
to study their poles and regions of convexity. 

We highlight in the following how the identification of reconstruction pairs (Definition l2.1D . and the 
analytical expressions for the Lagrange reconstructing polynomial (Proposition 14.51 * and its approxi- 
mation error (Proposition l4.7l l. can be used to develop general analytical expressions (valid VM± e Z : 
M :- M_ + M+ > 2 and VA's e |1, - - ,M - 1}) for the weight-functions, prove that there can be no poles 
at cell-interfaces (« + ^ Vn e Z), and extend the important results obtained in Liu et al. [19]. It is 
quite straightforward, using the definition of the Lagrange reconstructing polynomial (Definition l2.3D . 
to show by (6sA that the polynomial ah,M_,M^AO (|45g[l appearing in the representation ( I45dl ) of the La- 



grange reconstructing polynomial on the stencil S,;m_,m+ is the reconstruction pair, on a unit-spacing 
grid, of the corresponding polynomial a/.M-.M+.K^) < l45hl l appearing in the representation ( I45el l of the 
Lagrange interpolating polynomial on the same stencil 

a/,,M_,M+/(f) = [^(l;l)(a/,M_,M+/)] (^) a/.M^.M^.K^) = J^^ ah,M-,M^A^) dj] I y^^g^]^ ^ ' (^8) 

It is easy to prove by ( l68l ). using the mean value theorem for the definite integral 1I27L pp. 3 50-359], 



and the knowledge of the M roots of a/,M_,M+ AO ic^f.M-.M,An) = V« € {M_, • ■ • , M+) \ {^Ti^,!!) that 



a/7,M_,M+,f(n +2)^0 ^ y„ g ^ (69) 

and that all of the M roots of ah,M-.M+A.O are real. It can be shown that both {ahM-.M^AO^ ^ ^ ■ ■ ■ , M+}] 

anA.\afM_M^AO^ ^ ^ {-M_,--- ,M+}} form a basis of the (M+ l)-dimensional vector space of polynomials 
of degree < M in ^, Mm[^], and therefore, none of these polynomials is identically 0. We can work out 
several identities for the polynomials Q'/,,M_,M+,f(^), with corresponding identities for Q'/,M_,M+,f(^) because 
of ( l68l ). and show that an analytical expression for the rational weight-functions a-h^M_,M+,K„kM) given 
by the recurrenc 

Ac - 1 



minCATs-l./ts) 



4=max(0,/(:s-l) 

V/ts = o,- - < M- 1 (70) 



This analytical formulation, which only requires ( 45g[ l as input, is easily programmed in any symbolic 



computation package, and can generate the rational weight-functions VM± eZ : M - M^+M+ > 2. Since 



The recurrence relation <70t for K^>2 holds also for a"/,A/-,M+,s:s,fc(f)- 
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all of the M roots of any of the polynomials Q'/,,m_,m+,K^) are real, by ( iTOl ). we can show by induction that 
all of the poles of the rational weight-functions o-k^^.M^.K^kM) are real. These analytical results can 
then be used to extend the results of Liu et al. 1 19] VM e N>2, and for arbitrarily biased stencils in 
a homogeneous grid, providing at the same time simple symbolic computation routines for roots and 
poles. 

Going into further details and results is beyond the scope of the present work. We include however 
the following result. For the - [jj-level subdivision of the usual WENO stencils lfl9ll S, |^mj^_|^mj, 

we know from direct computation il9tl that the weight-functions o",, |^mj^_|^mj |-m-| ^ (i) > 0. Using the 
analytical expression ( l70l > we have obtained computationally the following result 

Result 6.1 (Positivity of linear weights at / + \). Assume that \M+\ < 9, satisfying ( l65] l. Then if for 
the subdivision level Ks o/'S,,m_,m+ ( 1651 * all substencils contain either point i or point i + 1 

then the rational weight-functions ( I70D satisfy 

(Th,M-,M^,K,,kXi) > e (0, ■ ■ ■ , K,] (71b) 

□ 



6.2. Truncation error o/WENO approximations to f'{x) 

Nonlinearity, ensuring monotonicity jssll, in WENO schemes is introduced by nonlinear weighting- 
out of stencils for which the reconstructing polynomial is nonsmooth. Smoothness is almost invari- 
ably IE i, [H El m, [H, III li, . 15] measured using the Jiang-Shu smoothness indicators LSJ. Let 



M : M — » K, M 6 N>i and Ax e R>o, and define 

M 



^M(x;Ax;M):=g J ^ ^ Ax^'-' \^—u(Oj d( (72a) 

By ( I72bl ) it is seen that, upon normalization by Ax, Pm{x\Ax;u) is [37] the usual norm of {d^uj in the 
Sobolev space '((-i, +j)) :- W'^ '-^((-5, +5)). The Jiang-Shu smoothness indicator, when consider- 
ing reconstruction at x,+ i on S;_m_,m+ (Corollarv l4.8l >. is define^H J^l, 4] by 

PpH,s,M-.M,,i+h Ax; phi-; S/,m_,m+, Ax)) (72c) 

with, as usual, M :- M_ + M+. The realization of the optimal order-of-accuracy by the WENO schemes 
hinges upon the fact fl, [H, EE Hi that j6p,„s,„ for different stencils of equal length M := M_ + 

M+ but different biasing around the pivot-point x, only differ to 0(Ax'^'^^), the lower-order part being 
common to all stencils of equal length M. Expansions in powers of Ax"/"/.^^" (n < [f J) have been 



The interval of integration was defined by Jiang and Shu 1 6] as the cell [-v,._^ , ^j^^], with center the pivot point .v; where we 
wish to approximate the derivative f- := /'(x,). This introduces some sort of upwinding, since the interval [.v,, could have 
been used (this would correspond to using /3m(.v^ u Ax; pi,(-; S,jvf__M+ ■ Ax)) instead in (72c). This is the choice made in the context 

of central WENO interpolation by Carlini et al. 136]. The choice of [.v,, .y,+ i] as the integration interval for upwind-biased WENO 
schemes has not been studied. Nonetheless, in the context of cell-centered finite-volume schemes |38] the choice of the cell as 
the volume of integration seems natural. 
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given in the literature, up to the WEN017 (composed by 9 substencils of length M - 8 cells using 
symbolic calculation la, 13, 14 , l^l- Using the analytical expressions for the error of the reconstructing 
polynomial with respect to the unknown function h(x) which is reconstructed, eg ( I56cl l. dSlcl * or dSldl l. 
it is quite straightforward to explain the existence of the common part. Using ( l56aD in ( l72cD yields, by 

dml ) 



4=1 ^-J 



(73) 



te the common (stencil-independent) part is indeed the Sobolev norm of d^hix; + ^Ax) and the non- 
common part involves the approximation error of the Lagrange reconstructing polynomial (Proposi- 
tion |4]7]l, and is therefore stencil-dependent. The deconvolution Lemma [2.51 is necessary to compute 
analytically the expansion of the common part in terms of powers Ax^f"f:^'" (n < [f J), and OSlcl l com- 
bined with the deconvolution Lemma [23] is used to evaluate analytically the expansion of the stencil- 
dependent part of ( [73] ). whose knowledge is essential for the evaluation and improvement [6, 13, 14] of 
the design of nonlinear weights, especially when interested in developing weights maintaining one of 
the great advantages of the Jiang-Shu weights, viz the straightforward extension to arbitrarily high- 
order accuracy 1 11, 15]. Furthermore these expressions were used to compute analytically the leading 
2 terms of the asymptotic expansions of the Jiang-Shu nonlinear weights [6] and from these the lead- 
ing term of the truncation error of WENO and WENOM tl3i] schemes. These developments are quite 
lengthy, and will be reported elsewhere. 

6.3. Extension to higher derivatives by multiple reconstruction 

The reconstruction approach can also be used to approximate f"(x) (and in general /*"'(x)), and in 
particular to compute interface fluxes, for high-order conservative discretization of diffusive terms in 
finite-volume methods, eg in the direct numerical simulation (DNS) of turbulence [39]. 

Assume the conditions of Definition l2.1[ for /(x) and h(x), but defined on / = [a- Ax, b + Ax] c M, and 
assume that 3 f) : / — > R which satisfies 



(Sa) 



1 r*x+ \ Ax 

x) = — " t)(Od( 

Ax Jj,_u 



Vx e [a - ^Ax, b + i Ax] 



(D 



h'-"\x) = 



+ iAx) - t)('"-"(x - iAx) Vx 6 [, 



A^ 



Vn 6 (1, 



J Ax, b + 5 Ax] 
■,N] 



(74a) 
(74b) 



assuming f,h, i) e C^[a - Ax,b + Ax], with e N>2. Recall (Remark 12. 8D that reconstruction pairs of 
continuous functions, if they exist, are unique. Combining ( [6a] ) with ( l74aD . we have 



/W = [^a;A.)(/^)](^) ^ [^aU) ° ^a;Ax) mix) = f^^' f T'^' 



drj 



(75) 



and we may write f) - /?(2;Ax)(/) the reconstruction pair of /(x) for the computation of the 2-derivative. 
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Differentiating ( [75l l twice with respect to x, we readily obtain the exact relation 



fix) = 



f)(x + Ax) - 2f)(x) + f)(x - Ax) 



1 

Ax 



Ax2 

f)(x + Ax) - \)(x) f)(x) - f)(x - Ax) 



Ax 



Ax 



(76a) 
(76b) 



•P/Ax+iAx) 'P/z'(x-iAx) 



By successive application of the deconvolution Lemma r2.5[ assuming the conditions of Lemma l2.51 we 
can obtain the deconvolution relations for ()(x) 



f=0 1.5=0 ^^■^ ^ ^1) 



^ U^n.itU^. + OfAx2L^J+2. Vx e [«, b] 



(77a) 



U=o 



AX / (x) + cy(Ax ) v„gi^,, .„<^_2L^J 



(77b) 



and redevelop all the results concerning R(i-ax) for /?(2;Aa), and in general for reconstruction procedures 
determining the interface fluxes for the computation of /"(x). In an analogous manner we can tackle 
the important practical problem of very-high-order conservative discretizations of {fAix)f^(x)y 1.4(^1 . 



7. Conclusions 

The results in this paper concern both the general relations between two functions constituting a 
reconstruction pair (Definition [211), and the analysis of the approximation error of the reconstructing 
polynomial (Definition I2.3D . 

We call a function h(x) whose sliding averages over a constant length Ax are equal to /(x) the 
reconstruction pair of /(x), h - R(i-Ax)(f) (Definition 12. ID . The exact relations Ax/'"^(x) - /z*"~''(x + iAx) - 
/;(n-i)(jj - I Ax) ^ are the basis of the numerical approximation of /'(x) by reconstruction f?,^,^^^. 

The reconstruction pair of the exponential function is [/?(i;Ax)(exp)](x) = gr(Ax)e '' (Theorem l2.9l >. The 
function gr(x) il9hi is the generating function of the numbers t„ (Tab. satisfying recurrence dlOcD . 
The numbers t„ ( I19cI > define the coefficients of the analytical solution (Lemma l2.5l > of the deconvolution 
problem for Taylor polynomials |2, (3.13), pp. 244-246]. This analytical solution (Lemma l2.5D is one of 
the main results of this work. It was also obtained by an alternative matrix-algebra oriented approach 
(B. 

The reconstruction pair of a polynomial of degree M e N is also a polynomial of degree M (Lemma l3.1l t. 
whose coefficients can be explicitly determined by ( I26fl ) using the numbers t„ (Tab. [T), R(i-^x) being a 
bijection of the vector space of polynomials of degree < M e N onto itself (Theorem lS.lD . 

The Lagrange reconstructing polynomial (Definition l2.3D on an arbitrary stencil S,.m_,m+ ■- {' - M^, ■ ■■ ,i + M+], 
on a homogeneous grid, in the neighbourhood of point i (Definition 14. Il l, is the reconstruction pair of 
the Lagrange interpolating polynomial on S,.m_.m+ '■- {i - A/-, • ■ • , / + M+). The Lagrange reconstructing 
polynomial on S,_m_,m+ is of degree M :- M_ + M+ (Proposition I4.5D and approximates h(x) to (9(Ax*^+') 
(Proposition 14.71 *. The complete expansion dSGaD of the approximation error of the Lagrange recon- 
structing polynomial in terms of powers of Ax can be expressed using the polynomials /l/,,M_,M+,«(f) 
defined by ( I56cl >. Most of the standard results of existence and uniqueness of the interpolating polyno- 
mial apply to the reconstructing polynomial (Theorem l5.3l *. 

Typical applications include the analytical expression and results on the roots and poles of the ra- 
tional weight-functions combining the Lagrange reconstructing polynomials on substencils of S,,m_,m+ ■- 
{i - M_, ■ • ■ , / + M+) into the Lagrange reconstructing polynomial on the entire stencil, the analytical 
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expression of the Taylor expansions of the Jiang-Shu |@] smoothness indicators and of the truncation 
error of WENO schemes, and the analysis of the discretization of ,f"\x) by «-reconstruction. It is hoped 
that the theoretical relations on reconstruction pairs and the analytical expressions of the approxima- 
tion error of the reconstructing polynomial will be useful in the analysis and improvement of practical 
discretization schemes. 

Appendix A. Useful relations for summation indices 



We summarize here several relations |22, 23] used in the text to manipulate the limits of summa- 



V« e Z 



+ 2k+l\ e / 2k+\\ 2k 



(A.1) 



tion indices, and some other useful formulas. 

a <n <=^ \a \ < n 

a <n <^=> \_a\ < n 
n</3 n < 1/31 

n</3 n < [J3\ 

s<2k <^ \^^<k 

s <2k <=^ Lf J < k 

2k <s <^ k< rf 1 

2k <s ^ k<lji 

^ ] ^ ^nm — ^ ^ ^n,s-n ~ ^ ^ ^ ^ ^ ^ ^s—m,m (A.3) 

1 l£ + 2k+l\ \ le + 2k\ 

(A.4) 



Ms^k&Z (A.2) 
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